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When climate change couples social neglect: malaria 
dynamics in Panama 

Lisbeth Amarilis Hurtado^, Lorenzo Caceres^, Luis Fernando Chaves^'^ and Jose E Calzada^ 

A major challenge of infectious disease elimination Is the need to Interrupt pathogen transmission across all vulnerable populations. 
Ethnic minorities are among the key vulnerable groups deserving special attention in disease elimination initiatives, especially because 
their lifestyle might be intrinsically linked to locations with high transmission risk. There has been a renewed interest in malaria 
elimination, which has ignited a quest to understand factors necessary for sustainable malaria elimination, highlighting the need for 
diverse approaches to address epidemiological heterogeneity across malaria transmission settings. An analysis of malaria incidence 
among the Guna Amerindians of Panama over the last 34 years showed that this ethnic minority was highly vulnerable to changes that 
were assumed to not impact malaria transmission. Epidemic outbreaks were linked with El Nino Southern Oscillations and were sensitive 
to political instability and policy changes that did not ensure adequate attention to the malaria control needs of the Gunas. Our results 
Illustrate how the neglect of minorities poses a threat to the sustainable control and eventual elimination of malaria in Central America 
and other areas where ethnic minorities do not share the benefits of malaria control strategies intended for dominant ethnic groups. 
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INTRODUCTION 

One hundred years ago, the Panama Canal was inaugurated, and ever 
since, it has been a major route for global commerce.' The Panama 
Canal was celebrated as a major achievement of civil and envir- 
onmental engineering at the time, and its construction was also a 
major achievement for public sanitation.^'^ Endemic yellow fever 
and malaria quickly waned after intense vector control measurements 
were taken in 'a hotbed of malaria'."* Malaria control in the Canal Zone 
(name of the Panama Canal under the colonial rule of the United 
States) emphasized the long-term benefits of environmental sanita- 
tion.*'^ This priority was derived from the need to render the Panama 
Canal suitable for the colonialist exploitation of a convenient trade 
passage.''' Nevertheless, wider health benefits were observed by all 
residents of the Canal Zone and its surrounding neighbor, the 
Repiiblica de Panama.^'^'^ Malaria was rare in the Panama Isthmus 
by the mid-1930s* and had almost disappeared by the mid-1950s.^ 

Parallel to the 'official' history of great achievements in malaria 
elimination from the Canal Zone and control in the Republica de 
Panama, there was a silent history of malaria transmission in 
Amerindian populations of the Panama Isthmus. Earlier epidemiolo- 
gical^'^ and entomological''" studies on malaria in the Canal Zone 
revealed that poor housing, human crowding and other consequences 
of social exclusion exacerbated the malaria burden among Amerindian 
populations. Nevertheless, following the initial assessment of the mal- 
aria situation in the Panama Isthmus, concerns about malaria in 
Amerindian minorities were largely absent in studies that claimed 
an almost final victory against malaria in Panama.^ 



Malaria transmission was considered to be negligible in the Panama 
Isthmus from the 1950s onwards, including at the time of the Torrijos- 
Carter treaty signature in 1979" — an agreement that warranted the 
return of the Canal Zone to Panamanian sovereignty'^ — and in 1999 
when the Panama Canal was returned to the Republica de Panama." 
In the 1950s, the Republica de Panama also began its journey toward 
fully recognizing the Amerindian ethnic groups as citizens, with the 
increased recognition of rights for the Amerindian ethnic groups 
through the creation of the San Bias (current Guna Yala) Comarca, 
the first of several Comarcas, i.e., territories designated to provide fuU 
autonomy and social integration for Amerindian minorities within the 
multicultural and ethnically diverse definition of the Republica de 
Panama."''"* In this context, Comarca Madungandi was created in 
1996. This region is part of the homeland of the Guna (formerly known 
as Kuna), an indigenous ethnic group of Panama (Figure lA). The 
Gunas are a Chibchan group. DNA analysis showed low variability in 
their mitochondrial sequences, which is compatible with a population 
bottleneck coincident with the time of ethnogenesis for the Chibchas, 
approximately 10000 years ago.'^ Madungandi is a 1200 km^ area 
within the Chepo district of the Panama Province, on the Atlantic basin 
in mid-eastern Panama. Puente Bayano, the most accessible village in 
Madungandi, is 80 km from Panama City and 120 km from Colon, the 
two major cities in Panama (Figure IB). In 2010, there were 4217 
Gunas living in the 14 villages that comprise Madungandi,*^ with most 
residing along the shores of the artificial lake Bayano (Figure IC). 
Approximately 90% of the Madungandi Guna natives live in extreme 
poverty, which is almost five-fold higher than the percentage observed 
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in the rest of Panama.'^ About half of the aduh population, i.e., 
people over 15 years of age, is illiterate (45.9%), and only 17.5% 
completed elementary education. The major economic activities are 
subsistence farming and fishing, which employ 56.4% and 17.2%, 
respectively, of the adult population.'^ The health status of the 
Guna lags behind that of other Panamanian citizens; the Guna have 
a life expectancy that is between 7 and 10 years shorter than that 
of other Panamanians, and a three-fold increase in the mortality 
rate for children under one year of age when compared with other 
Panamanians.'^ Malaria is the most common infectious disease 
affecting the Guna. Since 2000, approximately 45% of the malaria 
cases in the entire Repiiblica de Panama have occurred in Comarca 
Mudangandi."'"* The Guna are also commonly affected by respir- 
atory and digestive infections, which are likely exacerbated by a lack 
of access to basic services." For example, 76.3% of the population 
has no access to potable water, and 84.3% has no access to basic 
sanitation, e.g., disposal of wastewater. Housing conditions are 
destitute and insect-friendly (Figure ID), with most houses buUt 
of wood and with thatched roofs, eaves, and unscreened doors 
and windows that facilitate Anopheles spp. mosquito infestations.^ 



The disproportionate share of the malaria burden of the Repiiblica 
de Panama that is held by the Guna Amerindians of Comarca 
Madungandi not only highlights their socially outcast condition but 
also their likely enhanced vulnerability to malaria, a parasitic disease 
with transmission that is sensitive to climatic variabiUty.^""^^ In this 
study, we analyzed the dynamics of malaria transmission among the 
Gunas from Comarca Madungandi, using a 34-year monthly malaria 
case time-series from Comarca Madungandi. We assessed whether 
major malaria epidemics observed in the region in the early 2000s" 
represent a major dynamical regime shift, i.e., an abrupt change in the 
mean and/or variability of malaria incidence. We also assessed 
whether the malaria epidemics were related to temperature increases 
in the region or whether they were related to changes in climatic 
variability associated with large climatic events, such as the El Niiio 
Southern Oscillation (ENSO). We found that inter-annual cycles in 
malaria epidemics were associated with ENSO during 1980-2013 and 
that a dynamical regime shift observed at the end of 2002 might 
be related to policy changes in the structure of the malaria control 
program, in which personnel and funding for this program were re- 
assigned within Panama's Ministry of Health without adequate 
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Figure 1 The Gunas of Comarca Madungandi, Republica de Panama. (A) Guna children at Puente Bayano elementary school. (B) Comarca Madungandi, Republica 
de Panama. The map shows the location of Comarca Madungandi in relation to Panama City and Colon. (C) Guna Village at the shore of Lake Bayano. (D) Traditional 
Guna house with earthen floor, wood walls, thatched roofs and no anti-insect screening. The inset shows the typical eaves that facilitate Anopheles albimanus 
mosquito infestation. 
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preparation to tackle malaria transmission. We consider that full 
inclusion of the Gunas into Panamanian society is the major challenge 
for successful malaria control, and eventual malaria elimination, in the 
Repiiblica de Panama. 

MATERIALS AND METHODS 

Malaria and demographic data 

Monthly malaria cases for Comarca Madungandi, from January 1980 
to AprU 2013, were obtained from the East Panama Region Vector 
Control Department of Panama's Ministry of Health (Departamento 
de Control de Vectores de la Region Panama Este, Ministerio de Salud, 
Repiiblica de Panama). The time series mainly records malaria cases 
from a sentinel health post in Puente Bayano, the main village in 
Madungandi (Supplementary Figure SIA). The time series only 
included malaria cases with diagnoses that were confirmed by the 
examination of Giemsa stained blood smears that were prepared using 
the thick smear method* and examined under the microscope by 
trained staff at the Gorgas Institute or the Panama Este regional 
Hospital. All positive slides, and 10% of the negative slides, were con- 
firmed by the Public Health Central Reference Laboratory of the 
Gorgas Laboratory, following the national guidelines for malaria con- 
trol of the Ministry of Health of Panama. For the study period, both 
the sensitivity and specificity of malaria diagnosis were consistent and 
were close to 100%, highlighting the quality of the data. During the 
study period, more than 95% of the malaria cases were due to 
Plasmodium vivax, but during the 2002-2006 outbreak, approximately 
12% of the cases were due to Plasmodium falciparum}^ 

Annual population size estimates (1980-2013) for the Chepo dis- 
trict were obtained from the Direccion de Estadistica y Censo of 
Contraloria General de Repiiblica de Panama. Data from the Chepo 
district was rescaled with the assumption that Gunas in Comarca 
Madungandi steadily accounted for 10% of the total population, as 
they did in the 2010 census (Supplementary Figure SIB). 

Monthly malaria incidence per 1000 was estimated as the ratio 
between the number of malaria cases and population size estimates 
for each month, multiplied by 1000. 

Study site and climatic data 

The natural vegetation of Comarca Madungandi is a humid tropical 
rainforest, with a high and stable temperature of approximately 
26-27 °C and a unimodal rainfall pattern with a dry season (mid- 
December to mid-AprU) and a wet season (mid-April to mid- 
December). Rainfall, maximum temperature (TMAX) and minimum 
temperature (TMIN) time series from the nearest weather station to 
Madungandi were obtained from ETESA, the Panamanian National 
Electrical Company. Sea surface temperature 4 (SST4, also known as. 
El Niiio 4 index) time series, an index for ENSO, a global climatic 
phenomena associated with extreme weather in Central America was 
obtained from the United States National Oceanic and Atmospheric 
Administration Climate Prediction Center (http://www.cpc.ncep. 
noaa.gov/data/indices/ersst3b.nino.mth.81-10.ascri). The United States 
National Oceanic and Atmospheric Administration data were collected 
from the area delimited by 5TSf-5'S and 160'E-150°W of the Pacific 
Ocean. 

Statistical analysis 

Trends in the climatic variables. We began our analysis by determin- 
ing whether trends (increases in the average value) were present in the 
local climatic variables: rainfall, TMAX and TMIN. We employed 
locally weighted scatterplot smoothing, a non-parametric regression 



method, to fit trends in the climatic time series data. For robustness, 
we also employed empirical mode decomposition (EMD), a method 
in which a time series is decomposed by building oscillatory signals — 
intrinsic mode functions (IMF) — that are repeatedly subtracted from 
the time series and a residue, a nonlinear trend, can be obtained after 
the oscillatory components are removed.^^ Details about the EMD 
decomposition algorithm and a mathematically rigorous explanation 
are presented by Huang et al?^ 

Seasonality of malaria incidence. The seasonality of the malaria 
incidence time series was assessed using box diagrams for each 
month.^'^ 

Malaria incidence time series correlation structure and association 
with climatic variables. We explored the correlation structure in the 
malaria incidence time series by inspecting its autocorrelation func- 
tion (ACF), i.e., the time series correlation with itself through time, 
and the partial autocorrelation function (PACE), i.e., the correlation 
between consecutive time lags.^'' Information from the ACF and PACE 
was used to fit an autonomous nuU model, i.e., a model without 
covariates, of the malaria incidence data. The nuU model was then 
used to pre-whiten the climatic covariates SST4, Rain, TMAX and 
TMIN. Pre-whitening is a process in which any common autoregres- 
sive structure is removed (often referred to as filtering) from a time 
series to study patterns of association with a focal time series.^"* 
Residuals from the autonomous null model together with pre-whi- 
tened residuals of the climatic covariates were then used to estimate 
cross-correlation functions (CCF) of malaria incidence with each of 
the climatic covariates. We used the information from the CCFs to 
build a full non-autonomous model that considered up to two lags 
from each of the climatic covariates, selecting the lags with the highest 
correlation, based on the statistical significance of the association and 
the magnitude of the association. This full model was then simplified 
to avoid having more parameters than necessary to understand the 
dynamics of malaria incidence. We simplified the fuU model by a 
process of backward elimination, i.e., removing the least significant 
covariate in rounds.^* We used the Akaike Information Criterion 
(AIC) to choose the best model in each round of model simplification. 
AIC is a model selection criterion that can be used to choose the best 
model after a trade-off function between the number of parameters 
and likelihood is minimized.^'' Finally, in all cases, assumptions about 
model error were verified using standard procedures for time series 
analysis.^* 

Breakpoints and dynamical regime shifts. To explore the possibility 
that climatic variables with trends and malaria incidence in 
Madungandi had dynamical regime shifts, i.e., abrupt changes in their 
average and/or variance during the studied period, we employed a 
parametric test based on an F distribution that compares the nuU 
hypothesis of no change with that of a significant change. Briefly, 
the method by Bai and Perron^^ allows the simultaneous estimation 
of several breakpoints, i.e., time points when a dynamical regime shift 
occurs, via an iterative algorithm by which the number and temporal 
location of breakpoints can be estimated by minimizing the Bayes 
Information Criterion (BIG). BIG is a model selection criterion that 
guides the best model selection when a trade-off function considering 
the likelihood of a model, the number of parameters and data points 
employed in parameter estimation is minimized.^"* 

After the breakpoints were estimated, we fitted the full non- 
autonomous model to the malaria incidence regimes, i.e., time series 
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segments, which were defined by the estimated breakpoints. These 
models were simplified based on AIC minimization and assumptions 
of the error were checked using standard procedures?^ We then 
studied changes in the magnitude of climatic forcing in the malaria 
time series by comparing the time series model parameter estimates 
for the different segments of the split time series.^* 

Non-stationary patterns of association in the time-frequency 
domain. The time-frequency association of non-stationary time ser- 
ies, i.e., the association of cycles in time series over time, can be studied 
using continuous wavelet transforms.^^ Here, we used cross-wavelet 
coherence analysis to determine whether peaks at a particular fre- 
quency and time in SST4 correspond to similar peaks in malaria incid- 
ence and local climatic variables (Rain, TMAX and TMIN) of 
Comarca Madugandi, Repiiblica Panama. We also depicted the asso- 
ciation between malaria incidence and climatic factors by first decom- 
posing the time series into a trend and oscillatory components via an 



EMD. The power spectra of the resulting trend and IMFs were 
estimated via the construction of a periodogram, which is a graphical 
representation of the distribution of power (or variance) among dif- 
ferent frequencies, the inverse of the period of cycles per year.^"* A peak 
in the periodogram indicates a dominant frequency, whose inverse is 
the period of the dominant cycles in a time series.^" Given the period of 
the cycles in the trend and IMFs from the EMD, we estimated their 
CCFs with the climatic covariates, to untangle the specific associations 
between the global (SST4, El Nino 4 index) and local (rainfall, TMAX 
and TMIN) climatic covariates with inter-annual, seasonal and noisy 
cycles in malaria incidence. 

RESULTS 

Figure 2 shows the climatic time series for Comarca Madungandi. 
Figure 2A shows rainfall (R) during the study period. Figure 2B and 
2C shows maximum temperature (TMAX) and minimum temper- 
ature (TMIN), respectively. Minimum temperature has a clear trend 
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Figure 2 Time series data from Comarca Madungandi, Republica de Panama (January 1980-April 2013) and correlation functions (A) Rainfall. (B) Maximum 
temperature in Comarca. (C) Minimum temperature. Lines representthetrends obtained with LOWESSand EMD. (D) El Nino 4 index, also known as, SST4. (E) Malaria 
cases per 1000 people. (F) Malaria transmission seasonality. Boxes delimit ttie 25th to 75th data quantiles, and the lines inside the boxes show the median of the 
distribution for each month. (G) ACFof the malaria time series shown in panel A. (H) PACFofthe time series shown in A. CCF between residualsof a SAR (2, 12) model 
for malaria cases per 1000 people and pre-whitened: (!) SST4, (J) rainfall, (K) maximum temperature and (L) minimum temperature. In G-L, the dashed lines indicate 
the 95% CLfor the correlations that can be expected by random. CL, confidence limits; LOWESS, locally weighted scatterplot smoothing. 
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of close to a 3 °C increase from 1988 to 2013, a result robustly found 
with locally weighted scatterplot smoothing and EMD (Figure 2C). 
Figure 2D shows the SST4, an index of ENSO associated with 
Panama's climate. Figure 2E shows the malaria incidence time series, 
in which the largest epidemics occurred in 1987 and 2005. During the 
34 years of the time series, a total 4976 cases were reported in the 
Comarca Guna of Madungandi. Figure 2F shows the malaria season- 
ality. A higher incidence of malaria tended to occur during the dry 
season (December-April). Malaria incidence also shows a clear auto- 
regressive pattern, supported by a decreasing ACF (Figure 2G) con- 
firmed by significant PACF, coefficients at 1, 2 and 11 months of lag 
(Figure 2H). The ACF and PACF suggested that malaria incidence can 
be described by the following autonomous (i.e., without covariates) 
seasonal autoregressive (SAR) model: 

M, = /i + (;ii(M,_i-/() + 02(M,-2-/i) + 0i2(M-i2-/i) 

+ (j)l(l)l2{Mt^l3-^i)+(l>2<l>i2{M,_li-^l)+E, 

where M denotes malaria incidence, is the mean malaria incidence 
for different time lags and e, is an identical and independent normally 
distributed error. The maxima of CCFs between the residuals of the 
model in Equation (1) and the pre-whitened SST4 (Figure 21), R 
(Figure 2J), TMAX (Figure 2K) and TMIN (Figure 2L) suggest that 
malaria incidence was significantly associated with SST4 (4), R (0), R 
(5), TMAX (0), TMAX (17), TMIN (17) and TMIN (8), with the 
number inside the parentheses indicating the lag of the association 
in months. 

With the information from the CCFs, we proceeded to fit the 
following non-autonomous (i.e., with covariates), forced, SAR 
model: 

M, = ^t + (|)^{M,_i-^l) + (|>2{M,_2-^l■) 

+ 02012 (M,_ 14 - ^0 + aSST4,_4 + i3Rt (2) 
+ /?2R, - 5 + yiTMAX, + 72TMAX,_ 17 
+ ,5iTMIN,_8 + <52TMIN,_i7 + e, 



M, = ;t( + 01 (M,_ 1 - ^i) + 02 (M,-2 - /■') 
+ 012(M,_i2-/O+0i0i2(M,_i3-/O 

+ 02012 (M,_ 14 - aO + aSST4,_4 + /^Rt 
+ ,5iTMIN,_8 + <52TMIN,_ 17 + fit 

Parameter estimates for the model described in Equation (3) are 
presented in Table 1. 

The dynamical regime shift analysis showed that no breakpoints 
(and their associated dynamical regime shifts) occurred in the 
TMIN time series after a model such as the one described in 
Equation (1) was employed to find breakpoints using the Bai and 
Perron^^ method (Figure 3A). In contrast, a breakpoint was found 
when a similar procedure was applied to the malaria incidence time 
series, on December 2002 (95% confidence limits September 2001- 
February 2004), as shown by the minimization of BIC for one break- 
point in Figure 3B. There were major differences in the seasonality of 
malaria during the two periods separated by the dynamical regime 
shift. From 1980 to 2002, although epidemics were more common 
during December, January and February, no major differences were 
observed across the different months (Figure 3C). In contrast, from 
2003 to 2013, a clear seasonality was observed, with the number of 
cases peaking during February, and the number of cases was signifi- 
cantly higher at the end of the wet season in November, the dry season 
(December-March) and the beginning of the wet season in April 
(Figure 3D). 

To test whether there were significant changes in the association of 
malaria incidence with climatic variables in the regimes defined by the 
breakpoint, the model presented in Equation (3) was fitted to malaria 
incidence for the periods from 1980-2002 and 2003-2013. After a 
process of model selection, we found that malaria incidence for 
1980-2002 (Supplementary Table S2) was a second order autoregres- 
sive process and was also significantly associated with El Nino 4 index 
(SST4): 

M, = /i + 0i(M,_i-Ai)+02(M,-2-/t) + °'SST4f_4 + St (4) 



Hereafter, this model will be referred to as the fuU model. This Using model selection (Supplementary Table S3), we found that 
model was simplified using the AIC (Supplementary Table SI), and malaria incidence for the period from 2003-2013 was best described 
several terms were dropped resulting in the following final model: by a first order seasonal autoregressive process and was significantly 



Table 1 Parameter estimates for the best models of malaria incidence (per 1000 people) in Comarca Madugandi, Republica de Panama, with 
and without breakpoints 



Estimates (±s.e.) 



Parameter (Lai 


g-month) 


No breakpoint 


Breal<point (December 2002) 
Before After 


/< 


Mean (-) 


3.06±0.55* 


2.45±0.53* 


4.80±0.95* 


^1 


Autoregressive first order, AR (1) 


0.348±0.049* 


0.171±0.057* 


0.697±0.066* 


^2 


Autoregressive second order, AR (2) 


0.322±0.049* 


0.376±0.058* 




in 


Seasonai Autoregressive, SAR (12) 


-0.106±0.052* 




-0.248±0.087* 


s 


El Nino 4 (4) 


1.71±0.62* 


1.46±0.62* 


2.14±1.19 




Rainfall (0) 


-0.0043 ±0.0017* 




-0.008±0.003* 


Si 


IVllnimum temperature (8) 


0.475±0.179* 








IVllnimum temperature (17) 


0.579±0.182* 






6| 


Variance of the error 


15.61 


14.13 


16.16 



* Statistically significant (P<0.05). 
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Figure 3 Breakpoint analysis. BIC for breal<point selection in (A) minimum temperature and (B) malaria incidence (per 1000 people). The breakpoint in B 
corresponds to December 2002. Seasonal malaria incidence in (C) January 1980-December 2002 and (D) January 2003-April 2013. In C and D, the boxes delimit 
the 25th to 75th data quantiles, and the lines inside the boxes show the median of the distribution for each month. 



associated with SST4 and rainfall (P<0.05): 

M, = + 01 (M,_ 1 - + 'i^'n {M,-i2 - y) 
+ 0if 2 (M,_ 13 - aO + =<SST4,_4 + /«t + fi, 

Parameter estimates for Equations (4) and (5) are presented in 
Table 1. The split model described by Equations (4) and (5) outper- 
formed the model in Equation (3), with the former having an AIC of 
2148, whereas the AIC for the other model was 2159, with both models 
having a total of eight parameters. The results that are displayed in 
Table 1 confirm the observations from Figures 2C and 2D, where a clear 
seasonal structure was significant only for the 2003-2013 period. 
Assumptions of error normality, variance homogeneity and temporal 
independence were confirmed for the models presented in Table 1 and 
Equations (1) and (2). Table 1 also shows that mean incidence was 
higher during the 2003-2013 period (approximately 5 cases/1000 people, 
for one month) than in the 1980-2002 period (approximately 2.5 cases/ 
1000 people, for one month), regardless of the two regimes having a 
similar unexplained variance (approximately 15, which was also similar 
to the unexplained variance of the model without a breakpoint). The 
association between malaria incidence and El Nifio 4 (SST4) was positive 
and became stronger after the breakpoint, yet was not statistically sig- 
nificant for the 2003-2013 period, when a negative and statistically 
significant (P<0.05) association with rainfall was observed (Table 1). 



The differences between the significant covariates in the malaria 
incidence models with and without breakpoints (Table 1) may be 
related to the non-stationary association between climatic factors 
and malaria transmission. Non-stationary associations can be visua- 
lized through cross-wavelet coherency analysis. Therefore, we per- 
formed cross-wavelet coherency analyses between SST4 and malaria, 
rainfall, TMAX and TMIN. SST4 was significantly associated with 
inter-annual malaria cycles of four-year periods between 1980 and 
1995, and eight-year periods between 1995 and 2006 (Figure 4A). 
Around 1995, four-year cycles in SST4 and rainfall were associated 
(Figure 4B), and a similar pattern was also observed for TMAX, which 
was also associated with SST4 near 2005 (Figure 4C). Both TMAX and 
TMIN where seasonally (period of one year) associated with SST4 
(Figure 4C and 4D). 

To further tease apart the non-stationary association between mal- 
aria and climatic factors, we decomposed the malaria incidence time 
series using EMD. The signal was decomposed into a trend and two 
IMFs, and the power spectra of the signal components were estimated 
(Figure 5). The components from EMD were the trend (Figure 5A), 
the first IMF (Figure 5B) and the second IMF (Figure 5C). The trend 
(Figure 5A) had most of its power, i.e., variability, at inter-annual time 
scales, with clear peaks that correspond to cycles of several periods 
(i.e., 1/frequency of cycles per year), primarily 16 and 4 years, as well as 
shorter cycles of 18 months (Figure 5D). The first IMF (Figure 5B) 
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Figure 4 Cross-wavelet coherence analysis. Coherence between SST4 and (A) malaria cases per 1000 people from January 1980 to April 2013, (B) rainfall, (C) 
maximum temperature and (D) minimum temperature in Comarca Madugandi, Republica de Panama. In all analyses, a smoothing window of six months was used. 
The cross-wavelet coherence scale is from zero (blue) to one (red). Red regions in the plots indicate frequencies and timesforwhich the two series share variability (or 
power). The cone of influence (within which results are not influenced by the edges of the data) and the significant coherent time-frequency regions (P<0.05) are 
indicated by solid lines. 



had a power peak that corresponds to seasonal variabihty, or one cycle 
per year (Figure 5E). The second IMF (Figure 5C) had a similar 
amount of variability through all of the frequencies, indicating that 
the variability can be considered to be white noise because no fre- 
quency dominated the power spectrum (Figure 5F). CCFs of the trend 
from EMD with SST4 (Figure 6A), rainfall (Figure 6B), TMAX 
(Figure 6C) and TMIN (Figure 6D) showed that inter-annual cycles 
in malaria incidence were associated with El Niiio 4. By contrast, the 
seasonality of malaria incidence was not correlated with SST4 
(Figure 6E) or rainfall (Figure 6F), but was significantly correlated 
with TMAX (Figure 6G) and TMIN (Figure 6H); a similar pattern 
was observed for the white noise component, which was also not 
correlated with SST4 (Figure 61) and rainfall (Figure 6J), and was 
significantly correlated with TMAX (Figure 6K) and TMIN 
(Figure 6L). Thus, the combined outcomes of the wavelet analysis 
and the CCFs of the EMD decomposed malaria incidence time series 
indicate that inter-annual malaria cycles are associated with El Nifio 4 
and rainfall, whereas the seasonal dynamics are mainly associated with 
temperature changes. 

DISCUSSION 

The continued transmission of malaria among the Gunas in Panama 
highlights the challenges of tropical disease elimination in multicultural 
and ethnically diverse nations. Although malaria was recognized as a 
major problem in the Amerindian populations by baseline research 
conducted prior to malaria control in the Canal Zone,*'' little to no 
attention was directed to malaria in Amerindian populations of the 
Panama Isthmus once the disease was effectively controlled in the 
Canal Zone and major economic centers of the Republica de 
Panama.^'^'^'"^" Lake Bayano was artificially created on the Chepo river 
basin in 1975. At that time, potential risks for vector-borne disease 
emergence among the Gunas, a vulnerable community living in the area, 
were ignored.^' After the lake was created, the Gunas, who used to live 
near the river, settled on the shores next to Lake Bayano.'^ Lake shores 



have abundant sunlit pools, which are ideal larval habitats for Anopheles 
albimanus, the major malaria vector in most of Panama and Central 
America. '"'^''"^^ Indeed, shordy after the creation of Lake Bayano, the 
proliferation of 'water lettuce' (Pistia stratiotes) and its associated mos- 
quito fauna were linked with a major outbreak of Venezuelan equine 
encephalomyelitis in the area.^^ Thus, the emergence of malaria in 
Madungandi might be partially explained by the lack of a health impact 
assessment prior to the creation of Lake Bayano and the lack of prevent- 
ive vector management afterwards.''^^"^' Indeed, a similar pattern was 
previously observed in Panama along the Chagres River basin.* 
Following the creation of Lake Gatiin in the lower Chagres River, com- 
munities in the upper Chagres River were significantly less affected by 
malaria than those on the shores of Lake Gatiin during the 1930s.^'* 
Nevertheless, the malaria burden in Madungandi was most likely unre- 
corded until the 1980s, when the Puente Bayano Health post was estab- 
lished to serve the indigenous communities around Lake Bayano. 

At least three great malaria epidemics have struck the Gunas of 
Madungandi since 1980, including one that occurred from 1987 to 
1990. An El Niiio event occurred in 1987. In 1988, the use of 
dichloro-diphenyl-trichloroethane (DDT) as a regular tool for indoor 
residual spraying (IRS) in endemic areas was terminated.'"' This malaria 
control policy change coincided with a time of major social instability 
associated with the end of Manuel Noriega's rule and the US invasion of 
Panama,"*^ which greatly affected the overall fiinctionality of public ser- 
vices in the Republica de Panama.^^ During 1997 and 1998, a new 
malaria epidemic emerged, synchronous with major epidemics observed 
in other parts of the world, especially in East Africa^"'"*^^ and South 
America,''''^'' which, like Panama, were subject to one of the strongest 
El Nifio events on record.''^ The third large epidemic began in 2002- 
2003, which also coincided with an El Nino event.''* This last epidemic 
lasted until 2005-2006, mainly because of problems associated with the 
decentralization of Panama's Ministry of Health, where the malaria vec- 
tor control program suffered a shortage of funding and personnel neces- 
sary to manage focalized outbreaks quickly after their identification. 
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Figure 5 Malaria incidence (per 1000 people) EMD and spectral analysis of IMFs. (A) Residue; (B) first IMF; (C) second IMF. Periodograms for (D) EMD residue, 
which shows most of the powerat inter-annual time scales, i.e., less than one cycle per year; (E) first IMF, which has the highest power at a seasonal time scale, i.e., one 
cycle per year; and (F) second IMF, which has high variability across the spectra and therefore can be considered as white noise. In D-F, the blue lines are the 
confidence intervals for the power spectra. 



Our analysis supports an extended role for ENSO as a potential driver 
of inter-annual malaria epidemics in Madungandi, given the robust 
results derived from our different methods of analysis, which showed 
a close association between the inter-annual cycles of ENSO, measured 
by SST4, and malaria epidemics in Madungandi. One of the dominant 
cycle periods of the epidemics is approximately 4 (exactly 3.7) years, and 
its harmonics, i.e., cycles whose periods are a multiple of the dominant 
period,^" such as 16 years, are oscillation periods observed in ENSO.'''^ A 
remarkable 10-year study of malaria in eastern Panama revealed the 
occurrence of epidemics every 4 years in the Chagres river basin, with 
the largest epidemic observed in 1939,^'^ an El Niiio year.'''* 

Multidisciplinary studies in the Chagres river basin in the 1930s^'* 
showed that dry El Niiio years were associated with the dynamic gen- 
eration of larval habitats in Lake Gatiin shores and riverbeds along the 
Chagres River basin, ideal habitats for Anopheles albimanus, the main 
malaria vector around lakes in Panama, including Lake Bayano.^"'^'^" 
The negative association between rainfall and malaria incidence in our 
models suggest that a similar biological mechanism might be driving 
the recent malaria epidemics in Madungandi. Furthermore, our wave- 
let analysis indicates that rainfall patterns were associated with SST4 in 
Panama during the study period, supporting the hypothesis that 
ENSO may be linked to changes in malaria risk levels, via the dynamic 



generation of An. albimanus larval habitats by abnormal rainfall pat- 
terns. Mathematical models have suggested that dynamic larval hab- 
itat generation can lead to increased entomological vector-borne 
disease risk via the increase in mosquito density.^^ In addition, recent 
epidemics could have also been enhanced by the rising minimum 
temperatures in eastern Panama, which might be reducing the vari- 
ance in overall temperature that mosquitoes experience during their 
development, a condition known to maximize the fitness of exo- 
thermic organisms. This enhanced fitness pattern in stable environ- 
ments has been observed in non-anopheline mosquitoes,^^ and in 
malaria parasite development inside anopheline mosquitoes,^* thus 
providing a biological explanation for the associations we found 
between malaria incidence, temperature and SST4. 

Nevertheless, the increase in temperature likely played no role in the 
dynamical regime shift that was observed for 2002. First, temperature 
was not an important variable in the split model. Additionally, P. 
falciparum invaded Panama in 2002 after being absent from the 
Panama Isthmus since the 1970s.'^ Moreover, the circulating P. falci- 
parum parasites carried mutations that conferred resistance to chloro- 
quine, the main antimalarial drug used to treat the disease at that 
time."'^'*This invasion was likely mediated by the constant movement 
of the Gunas across their homeland, which comprises certain malaria 
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Figure 6 CCF between malaria incidence (per 1000 people) IMFsfrom an empirical mode decomposition, EM D and climatic factors. EMD residue and (A) El Nino 4 
index, also known as, SST4, (B) rainfall, (C) maximum temperature and (D) minimum temperature. First IMFand (E) SST4, (F) rainfall, (G) maximum temperature and 
(H) minimum temperature. Second IMFand (I) SST4, (J) rainfall, (K) maximum temperature and (L) minimum temperature. In all panels, the dashed lines indicate tlie 
95% CL for correlations that can be expected by random. CL, confidence limit. 



hotspots in northern Colombia where human mobility has been 
shown to play a major role in malaria persistence,^^ a possibility sup- 
ported by mathematical modeling.^^ An interesting characteristic of 
the shift was that mean incidence increased, but there was no change in 
the incidence variance. In contrast, dynamical regime shifts observed 
in East Africa in the late 1990s had breakpoints, which indicated 
changes in both the variance and mean incidence,^* which would be 
expected in the presence of a biological amplification of temperature 
trends by vectors. The control of the epidemic that began in 2002 
also illustrates the boldness necessary to contain malaria on the edge of 
elimination. In 2005, a Vector Control Task Group was created and 
funded by the Panama Ministry of Health to tackle the malaria epi- 
demic in Madungandi." This administratively autonomous group 
devised an intensive operational plan to guarantee the following activ- 
ities in the remote areas affected by malaria: active surveillance, rapid 
outbreak containment and a high IRS coverage, aiming to cover at 
least 80% of the domiciles in Madungandi.'^ Thus, by means of 
administrative modifications, without changes in the attention guide- 
lines or malaria control procedures, a significant drop in the incidence 
rate was observed 4 months after the installation of the program. At the 



end of 2005, the annual parasite index (API) was 1.4, which represents 
a 29% decrease compared with the API observed in the previous year 
(API=1.7%). This decreasing trend continued during 2006 with a 
70% reduction (API=0.5%) and during 2007 with a 76% reduction 
(API=0.4%). Additionally, the autochthonous transmission of P. fal- 
ciparum was eliminated and malaria ceased to be a cause of mortality 
in Madungandi. 

Some of the current challenges for malaria control and elimination 
among the Gunas are biological. Insecticide resistance has become 
common among An. albimanus in Madungandi, where this mosquito 
species is resistant to pyrethroids, but susceptible to organopho- 
sphates and carbamates,^* which has led to a change in the selection 
of insecticides, similar to when An. albimanus was among the first 
malaria vector species in which DDT resistance was observed,^'* with 
vector control now being based on the application of organopho- 
sphates and carbamates.^* In addition to An. albimanus, other vectors 
of concern in the area are An. punctimacula,^^ An. triannulatus and An. 
apicimacula, whose larval ecology is not related to lake shore habitats^* 
and whose control will not benefit from environmental modifications 
to the shores of Lake Bayano. However, the major obstacles to malaria 
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control and elimination among the Gunas are likely not biological, but 
socioeconomic.^" A lack of: (i) intercultural understanding, (ii) eco- 
nomic inclusion of the Gunas and (iii) multinational political com- 
mitment, heavily undermines any sustainable malaria control or 
elimination among this ethnic minority and, consequently, from the 
Repiiblica de Panama. Control measures such as IRS tend to be unwel- 
come among the Gunas and are probably not the most effective mal- 
aria control tool for their traditional housing style.'' Due to lack of 
consent, the average IRS coverage per village at best reaches 50%, with 
coverage ranging between 12% and 70%. Efforts for housing 
improvement require a better understanding of 'western housing' 
benefits by the Gunas, who are, in general, cautious of any policy that 
is supposed to improve their livelihoods," due to a long history of 
neglect and abuse by the dominant groups of the Panama Isthmus. 
Moreover, the deplorable socioeconomic exclusion status of the 
Gunas renders impossible housing improvement as an endogenous 
initiative, led and funded by an 'autonomous' Comarca. Adapting 
traditional Guna housing to be malaria proof requires multidisciplin- 
ary research, and excellent communication between researchers and 
the Guna communities to ensure a successful change, as shown for 
vector-borne diseases.*' Poverty and wariness of the government can 
also interact and lead to the misuse of other effective malaria control 
tools, as observed with insecticide treated nets in Lake Victoria 
Basin.*^ Similarly, the nature of the Gunas as a multinational ethnic 
group highlights the need for setting a common agenda for malaria 
elimination across borders, and focusing on the ethnic group beyond 
the nationality of its members. Similar policies have been suggested 
and implemented for mobile populations, including the Ngabe, 
another major Amerindian group in Panama, many of whom season- 
ally migrate across the Panama/Costa Rica border for economic rea- 
sons.^"* For the Gunas, the commitment requires the development of 
an integrated policy between Colombia and Panama, which is hin- 
dered by the ongoing civil conflict in the former.*' Thus, a major 
challenge for malaria control and elimination among the Gunas, 
and from the Republica de Panama, is the fuU inclusion of native 
Amerindian minorities into a multiethnic and pluricultural 
Panamanian society. 

In summary, our data showed a clear association between ENSO 
and the incidence of malaria in Comarca Madungandi, which put into 
context, is likely exacerbated by the precarious status of the Gunas 
within Panamanian society. Our analysis also supports that a dynam- 
ical regime shift in transmission that was observed at the end of 2002 
was likely related to malaria control policy changes, wherein the re- 
assignment of personnel and funding for malaria control within the 
Ministry of Health undermined the ability to combat malaria. This 
policy change also coincided with the re-introduction of P. falciparum 
in Panama, including drug-resistant parasites, and these elements led 
to an increase in malaria incidence among the Madungandi Gunas, 
which was controlled only after the malaria control unit was re- 
assembled. 
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